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Abstract 

First, we study geometric variants of the standard set cover motivated 
by assignment of directional antenna and shipping with deadhnes, providing 
the first known polynomial-time exact solutions. 

Next, we consider the following general capacitated set cover problem. 
There is given a set of elements with real weights and a family S of sets of 
elements. One can use a set if it is a subset of one of the sets on our lists 
and the sum of weights is at most one. The goal is to cover all the elements 
with the allowed sets. 

We show that any polynomial-time algorithm that approximates the un- 
capacitated version of the set cover problem with ratio r can be converted to 
an approximation algorithm for the capacitated version with ratio r + 1.357. 

In particular, the composition of these two results yields a polynomial- 
time approximation algorithm for the problem of covering a set of customers 
represented by a weighted n-point set with a minimum number of antennas 
of variable angular range and fixed capacity with ratio 2.357. 

Finally, we provide a PTAS for the dual problem where the number of 
sets (e.g., antennas) to use is fixed and the task is to minimize the maximum 
set load, in case the sets correspond to line intervals or arcs. 
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1 Introduction 



In this paper, we study special geometric set cover problems and capacitated set 
cover problems. 

In particular, the shapes of geometric sets we consider correspond to those of 
potential directional antenna ranges. Several geometric covering problems where 
a planar point set is to be covered with a minimum number of objects of a given 
shape have been studied in the literature, e.g., in OIHE]. 

On the other hand, a capacitated set cover problem can be seen as a general- 
ization of the classical bin packing problem (e.g., see [S]) to include several types 
of bins. Thus, we are given a set of elements {1, . . . , n}, each with a demand di, 
and a set of subsets of {1, . . . , n} (equivalently, types of bins), and the goal is to 
partition the elements into a minimum number of copies of the subsets (bins) so 
the total demand of elements assigned to each set copy does not exceed a fixed 
upper bound d. 

Capacitated set cover problems are useful abstraction in studying the problems 
of minimizing the number of directional antennas. The use of directional antennas 
in cellular and wireless communication networks steadily grows [H [T3l [T5l [H]. 
Although such antennas can only transmit along a narrow beam in a particular 
direction they have a number of advantages over the standard ones. Thus, they 
allow for an additional independent communication between the nodes in parallel 
[14j . they also attain higher throughput, lower interference, and better energy- 
efhciency P[l3l[l5]. 




Figure 1: The sectors correspond to the reaches of directional antennas. 

We consider the following problem of optimal placement of directional anten- 
nas in wireless networks. 

There is a base station coupled with a network infrastructure. The station 
transfers information to and from a number of customers within the range of 
directional antennas placed at this station. Each customer has fixed position 
and demand on the transmission capacity. The demands are unsplittable, thus a 
customer can be assigned only to a single antenna. One can choose the orientation 
and the angular range of an antenna. When the angular range is narrower an 
antenna can reach further so the area covered by any antenna is always the same. 
There is a common limit on the total bandwidth demand that can be assigned to 
an antenna. The objective is to minimize the number of antennas. 

Berman et al. termed this problem as MinAntVar and provided an approx- 
imation polynomial-time algorithm with ratio 3 p]. They also observed in [2] 
that even when the angular range of antennas is fixed, MinAntVar cannot be 
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approximated in polynomial time with ratio smaller than 1.5 by a straightforward 
reduction from PARTITION (see [8]). 

We provide a substantially better polynomial-time approximation algorithm 
for MinAntVar achieving the ratio of 2.357. Our algorithm is based on two new 
results which are of independent interest in their own rights. 

The first of these results states that a cover of the set of customers with the 
minimum number of antennas without the demand constraint can be found in 
polynomial time. Previously, only a polynomial-time approximation with ratio 2 
as well as an integrality gap with set cover ILP were established for this problem 
in [2]. 

The second result shows that generally, given an approximate solution with 
ratio r to an instance of (uncapacitated) set cover, one can find a solution to a 
corresponding instance of the capacitated set cover, where each set has the same 
capacity, within r + 1.357 of the optimum. 

Berman et al. considered also the following related problem which they termed 
as BinSchedule [2j. There is a number of items to be delivered. The i-th item 
has a weight di, arrival time ti and patience pi, which means that it has to be 
shipped at latest by ti + Pi- Given a capacity of a single shipment, the objective 
is minimize the number of shipments. 

Similarly as Berman et al. could adopt their approximation for MinAntVar 
to obtain an approximation with ratio 3 for BinSchedule [2], we can adopt 
our approximation for MinAntVar to obtain a polynomial-time approximation 
algorithm with ratio 2.357 for BinSchedule. 




Figure 2: The X coordinate of an item i encodes U and the Y coordinate encodes pi. 
Shipment has capacity 10. The numbers indicate the weights. Items which are to be 
shipped together must be enclosed by an angle. 

Our third main result is a PTAS for a dual problem to capacitated set cover 
where the number of sets (e.g., antennas) to use is fixed and the task is to minimize 
the maximum set load, in case the sets correspond to line intervals or arcs. In 
the application to directional antennas, the aforementioned correspondence comes 
from fixing the radius and hence also the angular range of the antennas and the 
problem has been termed as MinAntLoad in [2]. The task is to minimize the 
maximum load of an antenna. In ^ , there has been solely presented a polynomial- 
time approximation with ratio 1.5 for MinAntLoad. 

Organization: In Section [2] we present problem definitions and notations. In 
Section [3l we derive our polynomial-time dynamic programming method for the 
uncapacitated variant of MinAntVar. In Section HJ we show our general method 
of the approximate reduction of the capacitated vertex cover to the correspond- 
ing uncapacitated one. By combing it with the method of Section [3l we obtain 
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the 2.357 approximation for MinAntVar. Finally, in Section [5l we present the 
PTAS for MinAntLoad, or more generally, for minimizing the maximum load 
in capacitated set cover of bounded cardinality, in case the sets correspond to 
intervals or arcs. 

2 Preliminaries 

This section presents terminology and notation used throughout this paper. 
We use U to denote {1,2,..., n}. If Xj G M are defined for i G [/ and A C U, 

An instance of the set cover problem is given by family S of subsets of ?7 = 
{1, . . . , n}; a cover is C C S such that IJasC ^ = U. We minimize \C\. An instance 
of capacitated set cover also specifies di for i £ U; a capacitated cover is a family 
of sets C such that (i) for each A £ C there exists B £ S s.t. A C B, while 
d(A) < 1; (ii) Uyiec ^ = U. Again, we minimize \C\. 

If for each j £ U we define radial coordinates {rj ,0j), we define angle sector 

with radius bound as 

7^(r, a, S) ={j £U : rj < r and Oj = a + P with 0</3<6}. 

In MinAntVar as well as its uncapacitated variant, U is the set of customers 
with radial coordinates defined in respect to the position of the base station. 
This is a variant of capacitated (or uncapacitated) set cover where S consists of 
sets of customers that can be within range of a single antenna, i.e. of the form 
TZ{r,a, p{r)), where p{r) is the angular width of an antenna with radial reach r. 

The trade-off function p is decreasing; to simplify the proofs, we assume that 
p{r) = 1/r, we can change the r-coordinates to obtains exactly the same family 
of antenna sets as for arbitrary p. 

3 Uncapacitated cover by antenna sets 

To simplify proofs, we will ignore the fact that the radial coordinate has a "wrap- 
around". We also renumber the customers so 9i < 6i+i for 1 < i < n. Observe 
that if 9i = 9j and rj > rj then every antenna set that contains i also contains j, 
so we can remove j from the input. 

It suffices to consider only n(n -|- l)/2 different antenna sets. For such an 
antenna set A, let i = mm A, j = max A. If i = j, we denote A as = {i}, 

and ifi < j, we set r{i,j) = {6j — Oi)^^ and define j4[i,j] = Tl{r{i,j),9i, l/r{i,j)). 
(This definition is more complicated when the "wrap-around" is allowed.) Because 
A C j] we can use ^[i, j] in our set cover instead of A. 

We say that points i and j are compatible, denoted i'^j, if i < j and there 
exists an antenna set that contains If i = j then is obvious; if i < j 

then i'v'j = {i,j} Q A[i,j] = ri,rj < r(i,j). If i'^j, we define S[i,j] = {k : i < 
k<j}\A[iJ]. 

We solve our minimum cover problem by dynamic programming. Our recursive 
subproblem is specified by a compatible pair i,j and its objective is to compute 
the size of minimum cover C[i,j] of ^[i,^] with antenna sets. If we modify the 
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input by adding the points and n + 1 with coordinates (^i — 1, e) and {On + 1, e) 
then our original problem reduces to computing C[0,n + 1]. 

If S[i,j] = then Cij = 0. Otherwise, S[i,j] = {oq, . . . , am-i}, where 
Ofc < ttk+i for k = 0, . . . ,m — 2. 

We define a weighted graph Gij = {Vij, Eij,c), where Vij = {0, ...,m}, 
{k,£+ 1) £ Eij iff' ak'^ai and for an edge {k,i+ 1) we define the cost c{k,£+ 1) = 
1 + C[ak,a£]. 

Note that Gij is acyclic. Therefore, we can find a shortest (i.e., of minimum 
total cost) path from to m in time 0{\Eij\) = O(n^) [6]. Let d be the length of 
this path. We will argue that C[i,j] = d. 

First, we show a cover of ^[i,^] with d antenna sets. A path from to m in 
Gij is an increasing sequence, and a path edge (n, v) with cost c corresponds to 
a cover of {a^, a„+i, . . . , a^-i} with A[au, a^-i] and c — 1 antenna sets that cover 
S[au, ttv-i]- 

Conversely, given a cover C of S'[-i,j], we can obtain a path with cost \C\ in 
Gjj that connects with m. 

For 74[A:,£] G C, we say that ^ — A; is its width. To make a conversion from a 
cover C of <S'[i, j] to a path in Gj.j, we request that C has the minimum sum of 
widths among the minimum covers of S[i,j]. 

This property of C implies that if A[k, P\ then: 

k and l are not covered by C — {A[A;, £]} (otherwise we eliminate A[k, P\ from 
C or replace it with a set that has a smaller width). 

From this we can conclude that for each pair of sets ^], ^[/c', G C, where 
k < k', one of two following cases applies: 

1. £< k', i.e., A[k,i] precedes A[k',i']; 

2. f < i, i.e., A[k',i'] is nested in A[k,i]. 

Let 2? be the family of those sets in C that are not nested in others. Clearly 
T) can be ordered by the leftmost elements in the sets. Note that if 74[/c,^] G P 
then for some /, g, c, we have 

Of = k £ S[i,j], 
ag = i £ S[i,j], 

c — 1 sets of C are nested in ^[/c,£] and they cover S[i,j], 
{f,g + 1) is an edge in Gij with cost c, 
g + 1 = m or A[ag+i,£'] G V for some £'. 

These (/, g + 1) edges form a path that connects with m with cost \C\. 

Our dynamic programming algorithm solves the n(n + 1)/2 subproblems spec- 
ified by compatible pairs i,j in a non-decreasing order of the differences j — i. In 
the reduction of a subproblem to already solved subproblems the most expensive 
is the construction of the graph Gij and finding the shortest path in it, both take 
quadratic time. Hence, we obtain our main result in this section. 
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Q ^ 

for {U G W) 

while {U / 0) 

Q^0 

for (i E [/, with di non-decreasing) 
if {d{Q) + di<l) 

insert i to Q 

remove i from U and P 
insert Q to Q 

Figure 3: FFD, First Fit Decreasing algorithm for converting a cover into a ca- 
pacitated cover. 

Theorem 1 The uncapacitated version of the problem of minimum covering with 
antenna sets n points, i.e., the restriction o/ MinAntVar to the case where all 
point demands are zero, can be solved in time 0{n^) and space 0{n^). 

Previously, only a polynomial-time approximation algorithm with ratio two 
was known for the uncapacitated version of MinAntVar [2]. 

4 From set cover to capacitated set cover 

By the discussion in the previous section, it is sufficient to consider only O(n^) 
antenna sets in an instance of MinAntVar on n points. Hence, MinAntVar is 
a special case of minimum capacitated set cover. 

Since we can determine a minimum uncapacitated set cover of an instance of 
MinAntVar by ignoring the demands and running the dynamic programming 
method given in the previous section, we shall consider the following more general 
situation. 

We are given an instance of the general problem of minimum capacitated set 
cover and a minimum set cover of the corresponding instance of minimum set cover 
obtained by removing the demands. The objective is to find a good approximation 
of a minimum capacitated set cover of the input instance. 

4.1 Approximation ratio r + 1.692 

We obtain an approximation with ratio 2.692 for minimum capacitated set cover 
on the base of minimum uncapacitated set cover hi* by running a simple greedy 
FFD algorithm (see Fig. [3]). Our analysis of this algorithm in part resembles that 
of the first-fit heuristic for bin-packing [5l [7], but the underlying problems are 
different. 

Theorem 2 Given an instance of capacitated set cover on n elements and an 
approximation with ratio r for minimum set cover of the uncapacitated version 
of the instance obtained by removing the demands, a capacitated set cover of the 
input instance of size at most r + 1.692 times larger than the optimum can be 
determined in time 0{v?). 
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Proof. To analyze FFD, we introduce a "slack function" s{x), and we also apply 
it to elements that we cover using notation notation Sj = s{di). Slack function 
has the following two properties: 

® if d{Q) < 1 then s{Q) < 0.692; 

@ if while {U 7^ 0) loop produces i + 1 solution sets, say Qo, . . . ,Qe then 

Let Q* be the optimum solution. Property (D implies that we start with 
s{P) < 0.692|Q*1. Property ® implies that algorithm FFD produces at most 
\U*\ + d{P) + s{P) < 2.692|Q*| sets. It remains to prove ® and ®. 

We define intervals = {x : -^p^ < x < k = 1,2,..., and we use them to 
divide P into classes, P^ = {i E P : di E Ik}- Now we define the slack function: 

We also introduce r{x) = s{x)/x and = r(di); observe that 

s{Q) < d{Q) msiXi^Q ri; 
< < i for i G Pfc. 

To prove ®, we look for the maximum possible s{Q). lfd{Q) < 1 and s{Q) > ^, 
then for some oq € Q we have > 5 , hence oq G -Pi , so dag = ^ + e for some 
£ > and Sag = ^. 

It remains to find maximum possible s{Q^) where = Q — {ao}- Note that 
d{Q'^) = I - £, (thus QnPi = 0). If s(Q^) > g then for some ai G we have 
fai^^, hence ai € P2 and dai = ^ + e for some e > 0. 

We can repeat the reasoning with = Qi — {ai} and conclude that it contains 
02 G Pq, and then with = — {02} we can conclude that it contains 03 G 
-P42, etc. Subsequent terms contribute very little to the overall result, so we can 
approximate the maximum possible s{Q) as | + ^ + ^ + 42x43 ~ 0.69103. 

The proof of property ® is in Appendix A. 

Since \U* | < n, our simple algorithm can be implemented in time O(n^). □ 
4.2 Approximation ratio r + 1.423 

FFD algorithm achieves the worst case behavior if the sets of the optimum solu- 
tions have demands of the form ^+^5 • • •} and the uncapacitated 
cover U* has sets that either have very small d(U), or group together all elements 
with a particular weight. 

E.g., for U that contains elements with da = ^ + £, algorithm FFD creates 
one-element sets. We can improve the approximation by preceding FFD with a 
phase in which we attempt to create "better sets". 

If d{Q) < 1 and Q Ci Pi = 0, the maximum s{Q) is obtained with demands 
i + e, i + e, i + e, ^ + £,..., and this yields = l + l + ^ + + ■ ■ ■ 0.4231. 

We can achieve the same even if there exists a G Q fl Pi if we reduce Sa from 
^ by about 0.269, to about 0.231. Then we need to modify the algorithm so it 
produces sets with d{Q) + s{Q) > 1. This is not necessarily possible, after all, Q* 
may even contain singleton sets. For this reason, we add the third term to our 
amortization of sets. For a G Pi we define 

Qa is the set in Q* such that a E Qa, 
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Xa = l- d{Qa)] 

Va = 0.1905 -s{Qa- {a}). 

For a Pi we set Xa = 0. Clearly, x{P) + d{P) < \Q*\ while s(P) + y{P) < 
0.4222|Q*|. Thus it will suffice to produce sets such that {d + x + s + y){Q) > 1, 
and for that, we just need to modify the way we create sets that contain elements 
of Pi. 

Let us consider what we (nondeterministically) can do, and what we need to 
do. Consider a E Pi and assume that da — ^ = Xa = Ua = 0. Then we can find 
S £ S and Ac S such that a e S, Ac P - Pi and s{A) = 0.6905. However, it 
suffices to find A such that 5(^4) = 0.269, less than 40% of what we can do. 

If we increase da, Xa or by some 5, both what we can do and what we should 
do decrease by 5, hence the ratio decreases. 

We can find a good candidate for A by "guessing" S E 5 and running an ap- 
proximation algorithm for the knapsack problem [12] in which items are elements 
i & S — {a}, the weights are di, the values are Vi = di + Si. It suffices to have 80% 
approximation . 

When we find a set Aa that has the maximum value (as returned by the 
approximation algorithm), we form set Ba = Aa U {a}. We do the following 
"accounting trick". For each i £ A and 6 G Pi, if 6 ^ a and i £ Qb, then we 
increase Xb by \vi. Thus we achieve da + Sa + Xa + Va + \v{A) > 1, while for the 
remaining elements b G Pi the ratio of what "they can do" (maximum possible 
v{Ah)) to what "they need to do" (the difference 1 — db — Xb — yb) remains bounded 
by 40%. 

After creating Ba for each a € Pi we run FFD algorithm with the remaining 
elements. 

In this preliminary version we omit details how to implement this refined 
algorithm in time 0(n|5|). 

Theorem 3 Given an instance of capacitated set cover on n elements and an 
approximation with ratio r for minimum set cover of the uncapacitated version 
of the instance obtained by removing the demands, a capacitated set cover of the 
input instance of size at most r + 1.423 times larger than the optimum can be 
determined in polynomial time. 

4.3 Approximation ratio r + 1.357 

One can observe that algorithm FFD has worst performance when some peculiar 
combinations of demands occur in sets of the optimum solutions, in terms of our 
classes, the worst pattern is (Pi, P2, Pe, . . .). Our second algorithm has an initial 
phase that handles all sets with an element from Pi; we decrease the slack for 
elements of Pi and spend more effort forming the sets, so even with the smaller 
slack we can amortize the cost of each set of our solution. 

Intuitively, members of Pi were troublemakers and our added phase took care 
of that. 

Because knapsack problem has fully polynomial-time approximation schema 
we could run a version with, say, 99% approximation, and this would allow to 
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decrease the slack in Pi by almost 0.6903/2. This would give an approximation 
ratio of about 2 + 0, 7/2 = 2.35. However, at this point we get another worst case 
— with the pattern (^2, P2, P3, Pi, ■ ■ ■)■ 

We say that a € P2 is a troublemaker if for some Q we have a & Q & Q* 
and |Q n P2I = 2. Here both elements of Q n P2 are troublemakers, we call them 
siblings. 

Now we will describe how to add a second phase to the algorithm so that the 
case of sibling troublemakers will cease to be the worst one. At that point we 
will have two classes of worse cases: (Pi, . . .), because they are compatible only 
with approximation ratios that are at least 2.35, and (P2,P3, • • •). The worst of 
the latter is (P2, P3, P3, P7, . . .). One can see that the slack of the latter is almost 
like the slack of the worst case of FFD, except that we have replaced a demand 
from Pi with two from P3, ^ + £ with two | + £. Thus this slack is approximately 
0.6903 - 0.3333 = 0.357. 

The second phase is similar to the first phase: we "guess" a set S* € 5, elements 
a,b & S n P2 and we run an approximation algorithm to find B £ S — {a, b} such 
that d{B) < 1 — da — db, while we maximize s{B). For all possible guesses, we 
pick one with maximum d{B) + da + dh, form set the Q = B U {a,b}, insert Q to 
our solution and remove Q from P. We repeat it as long as there exists S E S 
with |5'nP2| > 2. 

After the second phase is completed, we finish by running FFD with the re- 
maining P, the set of still uncovered elements. 

To analyze the second phase wc introduce a negative slack for each pair of 
sibling troublemakers, 0.1. When we form a set that contains troublemakers, we 
amortize it with the sum of the demands and slacks of elements, plus the slacks 
(and extra terms) of the troublemaker sibling pairs that arc involved. 

One can sec that the sum of slacks in Q E Q* that has a pair a,b of trouble- 
makers is at most 0.323 — we specifically decreased it by 0.1. We also define the 
extra terms similarly as before: 

Xa,b = 5 - d{Q - {a,b}); 

ya,b = 0A23-s{Q). 

If Xa,b + ya,b = 0, then the pair a, b "needs to find" 0.1, and it "can find" 0.423, 
so it suffices if it finds 25% of what it can find. When Xa,b (or ya,b) is positive, 
it decreases the "need to find" and "can find" by the same amount, so the ratio 
only improves (decreases). 

Now suppose that we form a set, and in the competition of "guesses" the 
winners were some a, 6 G P2. The critical case is when they are both troublemak- 
ers, each with its sibling, a' and b' respectively, and needs, Na and Nb- Because 
a, a' could find 4:Na, b,b' could find 4:Nb, they could find at least the average, 
2(Na + Nb). By applying | approximation, they found at least |(-/Va + N + b), 
the use | of that to satisfy their needs, and | of that to compensate the trouble- 
makers whose now can find less. The compensated troublemakers maintain their 
25% ratio of need/can. 

In this way, we obtain our strongest approximation results. 

Theorem 4.1 Let an instance of capacitated set cover be specified by a universe 
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set P = {1, demands di > for each i ^ P, and a family S of subsets of 

P. If an approximation with ratio r for minimum set cover of the uncapacitated 
version of the instance (i.e., where the demands are removed) is given then a 
capacitated set cover of the input instance of size at most r + 1.357 times larger 
than the optimum can be determined in polynomial time. 

Corollary 4.2 There exists a polynomial-time approximation algorithm for the 
problem of MinAntVar with ratio 2.357. 

By the reduction of BinSchedule to MinAntVar given in [2], we also obtain 
the following corollary. 

Corollary 4.3 There exists a polynomial-time approximation algorithm for the 
problem o/ BinSchedule with ratio 2.357. 

5 PTAS for MinAntLoad 

In MinAntLoad problem, the radius of antennas is fixed and the number m of 
antennas that may be used is specified. The task is to minimize the maximum 
load of an antenna. In [2], there is presented a polynomial-time approximation 
with ratio 1.5. 

In the dual problem MinAnt, the maximum load is fixed and the task is 
to minimize the number of antennas. Recall that achieving an approximation 
ratio better than 1.5 for the latter problem requires solving the following problem 
equivalent to Partition. 

Suppose that all demands can be covered with a single set, the load threshold 
is D and the sum of all demands is to 2D. Decide whether or not two antennas 
are sufficient (which holds if and only if one can split the demands into two equal 
parts). 

However, in case of the corresponding instance of MinAntLoad, we can apply 
FPTAS for the SubsetSum problem [11] in order to obtain a good approximation 
for the minimization of the larger of the two loads. 

If all demands can be covered by a single antenna set (and the sum of demands 
is arbitrary) then MinAntLoad problem is equivalent to that of minimizing the 
makespan while scheduling jobs on m identical machines. Hochbaum and Shmoys 
showed a PTAS for this case in [lOj . 

Interestingly enough, the PTAS of Hochbaum and Shmoys can be modified 
for MinAntLoad, while it does not seem to be the case with their practical 
algorithms that have approximation ratios of 6/5 and 7/6 jlOj . 

Because radial coordinate does not matter in MinAntLoad, the input is a 
sequence of pairs {9i,di), i = l,...,n. Initially, we ignore the issue of "wrap- 
around" so the antenna sets are of the form TZ{a) = {j £ U : a < 6j < a -\- Q}. 
Without loss of generality we assume that U = {1, . . . ,n} and 9i < 62 < . . . < On- 

In our PTAS, we try different values of the maximum load D. We can start 
using simple factor 2 approximation and then we can apply binary search. We 
will find an exact solution for a transformed problem in such a way that (a) the 
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cost of the optimum cannot increase, (b) a solution for the transformed problem 
can be converted to an actual solution while increasing the cost by a factor of 
l + e. 

For a fixed k, we will describe an (1 + e)-approximation algorithm that runs 
in time 0(n'^+'^), where c is a universal constant, while e ^ {1 + \nk)/k. 

We start by defining thresholds tj = D{1 + eo)^* and classes: 
Ci = {j £ U : tj+i < di < ti}, i = 0, . . . , k — 1 (large demands) and 
Ck = {j & U : di < tk} (small demands). We also set ei = and e = £o + ei- 
One can show that e is minimized when ei ^ 1/k and Eq ^ Ink/k. 

We will find exact solution to a problem where we have the same input but we 
re-define the cost/load of sets so (a) it cannot decrease and (b) if the new cost of 
Q satisfies cost{Q) < D then d{Q) < {l + e)D. We caU this problem Decreased. 

Intuitively, we divide elements into small and large. In the case of large el- 
ements, with dj > t, we decrease dj to d'j to have a small number of distinct 
values. In the case of small elements, we want to apply "greedy packing" and we 
"decrease" their contribution by not counting the last of them. More formally, we 
define decreased/relaxed instance Decreased as follows: 

for j £ Ci, we set d'^ to U+i, 

if Q n Ck = 0, we set cost{Q) to d'{Q), i.e., X^jgg d'^, otherwise 
if j = max((5 n Cfc), we set cost{Q) to d'{Q — Ck) + d{Q fl Ck) — dj, 
the task is to minimize maxg^g cost{Q). 

Clearly, the optimum of our Decreased instance cannot be larger than the 
optimum for the initial MinAntLoad instance. Also, since if j € Cj for i < k 
then dj < {l + €o)d'j and otherwise dj < d'j + e we conclude that d'{Q) < D implies 
d{Q) < {l+e)D. Thus, an exact polynomial-time algorithm for Decreased yields 
a PTAS for MinAntLoad. 

We say that a partition Q oi U is ordered if we have the following implication: if 
Q, Q' G Q, max(Q) < max(g'), j e Q n d, j' e Q' n d, then j < /. 

Lemma 1 For every solution Q' of MinAntLoad there exists an ordered solu- 
tion Q o/ Decreased such that maxQ^QCost{Q) < ioislkqqq' d{Q). 

Proof. We can transform Q' to an ordered Q in such a way that during that 
process for every Q £ Q' we will preserve \Q (1 Ci\ for each i > k and we will 
not increase d{Q n Cfc). Before Q is "finalized" we will allow fractional values for 
statements [j G Q] if j e Ck- 

Consider Q Q' that has minimal max{Q) and suppose that there exists 
Q' e Q' - {Q} and € d, j < f such that [j € Q'] > and [j' € Q] > 0. 
If i < k, we move j to Q and / to Q'; this does not change cost{Q — Ck) and 
cost{Q' — Ck). li i = k, let x = min{[j G Q'], [j' G Q]}, we increase [7 G Q] and 
[j' G Q'] by X and we decrease [j' G Q] and [j G Q'] by the same amount. This 
does not change d{Q n Ck) and d{Q' nCk)- 

When such Q', i,j,j' do not exists, suppose that there exists j G Ck such that 
< [j G Q] < 1; in this case j = max(Q fl Cfe); we increase [7 G Q] to 1 and for 
Q' ^ Q we decrease [j G Q'] to 0. This does not increase cost{Q) because cost 
does not count the last small element in Q. 
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Now Q and any other Q' satisfy the condition of ordered and we can remove Q 
and its elements from further consideration — and insert Q to Q. We repeat this 
until all sets are removed from Q' . □ 



The algorithm based on the lemma can be as follows. We represent a partial 
solution as counts (cq, . . . ,Cfe), that mean Cj elements of class Cj were covered. 
The are at most n^^Q|Cj| < {n/{k + 1))'^"''^ such partial solutions. Because we 
add sets to a solution in order of increasing max(Q), a partial solution covers Cj 
smallest elements of Cj — smallest in terms of their j's, or, equivalently, Oj^s. 

Adding a set to a partial solution (cq, . . . , c^) is an edge to another such vector, 
(cq, . . . , c'^). Such an edge is determined by the sequence (cq, . . . , c'^_i), because 
then we can find maximum possible . An edge is valid if it implies the increase 
in the maximum index of a covered element, and Y^\Zo{c'i — Ci)ti+i < D. Because 
a new set can cover at most l/t^ = k large demands, the number of possible edges 
is below 4^. We need to find the shortest path from (0, . . . , 0) to (|Co|, . . . , \Ck\), 
and we can use breadth first search; thus the time is proportional to the number 
of edges, or 0((4n//c)''+^). By e ^ {l + lnk)/k, the time can be also expressed as 
,^iin7+0(i)^ Hence, we obtain our PTAS for MinAntLoad. 

Theorem 4 MinAntLoad forn points admits an approximation with ratio 1+e 
in time n~^^~~^^^^\ 

Note that the only geometric property of antennas with fixed radius that we 
used to design the PTAS for MinAntLoad is their correspondence to intervals 
or arcs. Hence, we obtain the following generalization of Theorem HI 

Theorem 5 The problem of minimizing the maximum load in a capacitated set 
cover where the sets correspond to intervals or arcs admits a PTAS. 

6 Concluding Remarks 

We are quite convinced that our general method of approximating with ratio 
r + 1.357 minimum capacitated set cover on the base of an approximate solution 
with ratio r to the corresponding minimum (uncapacitated) set cover can ulti- 
mately achieve the ratio r + 1.3. In particular, this would improve the ratio for 
MinAntVar to 2.3. It seems however that some new ideas are needed to obtain, 
if possible, ratios below r + 1.3 and 2.3, respectively. 

Our aforementioned method can be also used to approximate optimal solu- 
tions to the natural extension of MinAntVar to include several base stations by 
combining it with known approximation algorithms for geometric set cover in the 
plane, e.g., B^M- 
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APPENDIX A: proof of the property © 



(D if while (U 7^ 0) loop produces i + 1 solution sets, say Qo, . . . , then 

We prove @ as follows. We remove from consideration every set Q created 
during that loop if d{Q) + s{Q) > 1. For j < £ we can define positive deficit 
Sj = l- d{Qj) - s{Qj). 

The claim is trivial if £ = 0, i.e. the loop creates only one set. Moreover, 
d{Qi-i) + d{Qe) > 1, hence it suffices to show s{Qi-i) + s{Q£) + Yl^^i^iQj) + 
s{Qj)) > £-1, equivalently, s{Qi_i) + s{Qi) > Yl^Zl Sj- 

Let tj be the time when algorithm FFD initializes ^ 0; and let -fL(j) be 
the class of the largest element of U time tj. 

If |C/ n PL{j)\ ^ -^(j) time tj, the algorithm would insert L(j) elements of 
PlU) to Qj, as each a G Pj^q) satisfies da + Sa > l^jj+t + L(i)(L(j)+i) = llj)^ ^^is 
would lead to in d{Qj) + s{Qj) > 1; a contradiction because we removed such sets 
from consideration. Hence \U nPii^j-^] < L{j) at time tj and the algorithm inserts 
entire remaining -Pl^) to Qj as well as at least one smaller element. This shows 
that L{j) is increasing with j. 

We will estimate the size of deficits and the "surplus" + s{Q(). 

First, we estimate s{Qj) in terms of A = L(j ' + 1). While we form set Qj, we 
can always insert an clement from P^, unless 1 — d{Q) < so Qj has a subset Q' 
with d{Q') > 1 — = and min^gQ' r,,, > j^, hence s{Qj) > x{x+i) — Est{X). 
Est{X) is decreasing with A, starting with A = 3. The case of A < 2 is not possible, 
because it implies that Qj has an element of Pi, hence, 6j < 0. 

Second, we apply the same reasoning for Qe-i U Qe and A = P{i): at time t^ 
there exists b € Pa and Qe-i U {b} contains a subset Q' such that d{Q') > 1 and 
miuaeQ' ra > j^, hence s{Qi^i) + s{Qi) > s{Q') > j^. 

Third, because we could insert b when we were creating Qj for j < i we have 
d{Qj) > 1 - i- 

Fourth, for > 1 we estimate S^^^; because A = L(i — A; + 1) < A — A; + 1, we 
have s{Qi^k) > Est{A - k + 1) = (A-fc+i)7A-fc+2) ' ^^'^^^ 

A, , <^ - (^ - 1) A-fc - (A-fc)(3-fc)+2 _ ^ 

oe-k ^ i a) (A_fc+i)(A_fc+2) ~ A(A-A;+l)(A-fc+2) —^^'^K'^)- 

Because A > 3, A — > 2, this shows that we have positive deficits only for 
k = 2,3 (for k = 1 the estimate refers to Qe-i and this set contributes to the 
surplus). Thus it suffices to show that — est{2) — est{3) > 0: 

1 A— 2+2 2 11 2 

A+T ~ A(A-1)A ~ A(A-2)(A-1) ^ A+T ~ (A-l)A ~ (A-2)(A-1)A ^ 
_J_ 1 

A+1 (A-l)(A-2) 

In our fourth point of the reasoning we observed that A — > 2, and the smallest 
value of fc is 2, so A > 4 and the above estimate is indeed positive. 
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